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ABSTRACT 

We study the reliability of dark matter (DM) profile determinations for rich clus- 
ters of galaxies based on weak gravitational lensing measurements. We assume that 
cluster DM density profiles are fully described by two parameters: a scale radius r s 
and a dimensionless central density 5 C . We consider three functional forms for the 
mass distribution: a softened isothermal sphere (SIS), and two models proposed on 
the basis of numerical simulations of CDM universes: the Navarro, Frenk, & White 
model (NFW) and the Klypin, Kravstov, Bullock, & Primack model (KKBP). A max- 
imum likelihood parametric method is developed and applied to simulated images of 
background galaxies at z—2.0 lensed by clusters at z — 0.2, 0.5, and 0.8. The images 
have field size (~ 10 arcmin) and angular resolution (~ 0.5 arcsec) characteristic of 
the new generation of 8-10 meter-class ground-based telescopes. The method recovers 
the correct virial mass of the DM halo, M200, to within ~ 10% at all redshifts. The re- 
covered values of S c and r s are only accurate to ^25%; however, their correlated errors 
are such that greater accuracy is obtained for the mass. Clusters with an SIS density 
profile are recognized by the method, while NFW and KKBP clusters are recognized 
as not being SIS. The method does not discriminate between NFW and KKBP pro- 
files for clusters of a given mass, because they are too similar in their outer regions. 
Strongly lensed arcs near the cluster centers would be needed to distinguish NFW 
from KKBP. The accuracy with which S c and r s are determined, and the different 
models are distinguished, improves with increasing redshift because the gravitational 
potential is sampled out to a larger radius for a given field size. These preliminary re- 
sults conclude that our maximum likelihood method is a promising candidate for both 
virial mass measurements and testing numerical predictions of DM halo structures 
although further tests that include real world effects still remain to be done. 



1 INTRODUCTION 



The study of rich clusters of galaxies, over a range of red- 
shifts out to z ~ 1, has emerged as an important tool for 
studying structure formation and constraining cosmological 
parameters (White, Efstathiou, & Frenk 1993; Bahcall & 
Cen 1993; Carlberg et al. 1997; Bahcall, Fan, & Cen 1997; 
Riechart et al. 1998; Eke et al. 1998; Girardi et al. 1998; 
Borgani et al. 1999; Blanchard et al. 2000). The theoreti- 
cal underpinning of this approach is the Press & Schechter 
(1974; PS) approximation for the abundance of collapsed 
objects as a function of mass and redshift in a given cosmo- 
logical model. The PS formalism is extremely sensitive to 
cluster virial mass since the fraction of collapsed objects of 
mass > M is given by the integral of the tail of a Gaussian 
function with the lower limit of integration dependent on 



M and cosmological parameters. Therefore, even relatively 
modest errors in virial mass measurement can drastically 
affect the derived cosmological parameters (Willick 2000). 
Accurate virial mass determination is thus of paramount 
importance. 

Observational data for clusters do not, however, directly 
yield their masses. Instead, they measure related properties 
such as X-ray luminosity and/or temperature, cluster rich- 
ness, or velocity dispersion. To date, efforts to constrain cos- 
mological parameters using clusters have been based mainly 
on either X-ray temperature/luminosity (Girardi et al. 1998; 
Donahue et al. 1998; Sadat, Blanchard, & Oukbir 1998; 
Riechart et al. 1998; Borgani et al. 1999; Blanchard et al. 
1999), or velocity dispersion (Carlberg et al. 1998; Tran et 
al. 1999) measurements as mass indicators. However, the 
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accuracy of such mass estimates is open to question. They 
depend on assumptions of hydrostatic equilibrium of the hot 
intracluster gas and virial equilibrium of cluster galaxies, re- 
spectively, which may not hold in a given cluster. It is thus 
important that alternative mass estimators, independent of 
the assumption of hydrostatic or virial equilibrium, also be 
applied when possible. 

Weak gravitational lensing of background sources by 
clusters is perhaps the most attractive such alternative (see 
Narayan & Bartelmann 1999, hereafter NB99, and Bartel- 
mann & Schneider 1999, hereafter BS99, for recent reviews). 
While strong lensing, in which large arcs or multiple images 
of a single source are produced, can provide precise mea- 
surements of the mass in the inner regions (r ^ 200/i -1 kpc) 
of clusters (e.g., Tyson, Kochanski, & DelP Antonio 1998), 
it does not constrain the virial mass of clusters, which typi- 
cally corresponds to radii of \~2h~ 1 Mpc. In contrast, weak 
lensing, which consists of small but coherent distortions of 
background galaxy shapes, can be detected statistically out 
to several h~ x Mpc. Cluster mass estimates have now been 
obtained from weak lensing data by a number of groups 
(e.g., Tyson, Valdes, & Wenk 1990; Kaiser & Squires 1993, 
1996; Smail et al. 1997; Squires et al. 1997; Hoekstra et al. 
1998; Luppino & Kaiser 1997). These estimates have gen- 
erally been found to agree well with the X-ray and velocity 
dispersion mass estimates (Lewis et al. 1999; Willick 2000). 

Several of the these weak lensing studies were done with 
the Hubble Space Telescope (HST) whose superb imaging 
capabilities ensure that it will play an important role in fu- 
ture projects as well. However, the primary observational 
driver for future work is likely to be improvements in the 
capabilities of the new generation of large aperture ground- 
based telescopes. Wide field (~ 10'-1°) imaging CCD cam- 
eras have recently or will soon come on line (see, e.g., Lup- 
pino, Tonry, & Stubbs 1998; Kaiser, Tonry, & Luppino 2000) 
at telescopes with excellent image quality (,$,0.5" seeing), 
such as CFHT, Subaru, Magellan, VLT, and Gemini. These 
instruments are extremely well-suited to weak-lensing stud- 
ies. In particular, they will fully cover the field of moderate- 
redshift (0.2 £ z £ 1) clusters, with data that probe out to 
and beyond the virial radius. 

The problem of cluster mass "inversion" , i.e. reconstruc- 
tion of the mass as a function of position from weak lens- 
ing data, has received considerable attention in recent years 
(e.g., Kaiser & Squires 1993; Seitz & Schneider 1997; Bartel- 
mann 1998). Most of these efforts have aimed at nonpara- 
metric mass reconstructions, in which no a priori assump- 
tions are made about the underlying mass distribution. How- 
ever, there is evidence from numerical simulations that the 
dark matter distribution in virialized haloes might be well 
described by simple parametric models. Navarro, Frenk, & 
White (1997; NFW) have proposed a model (§2) that has 
received perhaps the most attention. Recent high-resolution 
simulations have suggested departures from the NFW model 
at small radii, leading Kravstov, Klypin, Bullock, & Pri- 
mack (1998; KKBP) to propose an alternative parameteri- 
zation. Recent efforts to distinguish between these alterna- 
tives using the rotation curves of dwarf galaxies have been 
inconclusive (van den Bosch & Swaters, 2000), leaving it an 



open problem. The NFW and KKBP distributions are, how- 
ever, similar at radii encompassing £ 10% of the virial mass, 
where both may be described as locally power-law profiles 
p(r) oc r~ 7 , where 7 varies slowly from $2 near a charac- 
teristic radius r s to 3 at r > r s . In this sense, NFW and 
KKBP differ distinctly from the generic softened isothermal 
sphere (SIS) model, in which p(r) oc r~ 2 at all radii beyond 
a small inner core. The SIS model yields the flat rotation 
curves characteristic of spiral galaxies, and thus has often 
been viewed as a reasonable representation of equilibrium 
dark halos. 

Given the findings by NFW and KKBP that dissipa- 
tionless gravitational clustering may produce simple and 
universal mass distribution laws, a strong case can be made 
for using weak lensing to constrain the parameters of such 
laws. This approach would complement the existing methods 
of nonparametric mass inversion. Fitting parametric mod- 
els has advantages and disadvantages. The chief advantage 
is that, if the model is correct or nearly so, a parametric 
approach is robust against random error in the presence of 
noise. The obvious disadvantage is that the method may not 
flag an incorrect model. Nonetheless, we believe that the 
case for universal density profiles is strong enough that it 
amounts to a "theoretical prediction" which can and should 
be fit to the data. 

This paper has two objectives, to see if we can distin- 
guish between several theoretically motivated density pro- 
files using weak lensing data; and second, to ascertain how 
accurately we can measure virial masses from such data. We 
adopt three trial parameterizations of cluster masses: SIS, 
NFW, and KKBP (see § 2 for details), and simulate lensed 
fields for each of these mass distributions, for clusters of 
varying mass and redshift. We then apply our maximum 
likelihood method to these simulations (§ 3). For each clus- 
ter of a given profile type, we compute the likelihood with 
respect to the parameters of that profile type and of the 
other two, allowing us not only to constrain the parameters 
but to test whether we can discriminate among the models. 
One of our key conclusions is that we can indeed distinguish 
the SIS from NFW or KKBP, but not NFW from KKBP. 
This is important, because a basic prediction of the Cold 
Dark Matter (CDM) hypothesis is that halo mass distribu- 
tions depart from an r~ 2 law in their outer regions. Our 
work shows that this prediction is testable via weak lensing 
data. 

In addition to developing a method to test theoretical 
predictions of cluster mass profiles, our work has another, 
more basic, goal, already alluded to in the first paragraph. 
For clusters to be used as cosmological probes via the PS or 
related formalisms, one must determine their virial masses 
with high accuracy. In the PS formalism, the term virial 
mass has a very specific meaning: it is the mass contained 
within a sphere whose mean interior density is a particular 
multiple, Sv, of the mean density of the universe. The quan- 
tity 8v depends on cosmological parameters and on redshift. 
Since the former are not known a priori — and indeed may 
be the goal of our study — one needs to be able to express 
mass as a function of radius in a flexible way. This cannot 
be done independently of a model for the density profile of 
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the cluster; moreover, once a model is adopted, its denning 
parameters must be reasonably well-known to obtain a virial 
mass for a given choice of cosmological parameters (see, e.g., 
Willick 2000 for a detailed discussion). 

The outline of this paper is as follows: Section 2 and 
3 introduce the basic lensing formalism and our paramet- 
ric maximum likelihood inversion method. We then discuss 
the simulations used to test the method in Section 4, and 
Section 5 presents the results of these simulations. Section 6 
discusses the basic conclusions and directions for further re- 
search. A derivation of the shear and magnification for SIS, 
NFW and KKBP clusters is in Appendix 1, whereas Ap- 
pendix 2 considers the computation of the virial mass for a 
cluster with known structural parameters. 



2 BASIC FORMALISM 

In this section we derive expressions for lensing angle, con- 
vergence, and shear for two-parameter, spherically symmet- 
ric mass distributions such as SIS, NFW, and KKBP. Such 
mass profiles are characterized by a density scale p c and a 
scale length r s * We follow NFW in writing p c — p cr itS c , 
where p cr it = 3Hq /87rG, and thus write the density law as 



p(r) = p crit 5 C r)(r/r s ) , 



(1) 



where the function r/(x) characterizes the shape of the den- 
sity profile. For the three mass distributions considered in 
this paper we have 



rj{x) 



(1 + * 2 )- 1 

x-°- 2 (i + x 2 y 



(SIS) 
(NFW) 
lA (KKBP). 



(2) 



Figure 1 shows these three profiles for the case of a clus- 
ter with virial mass M200 = 10 15 /i -1 Mq and virial radius 
^200 = 1.63/i~ Mpc, for S v — 200. Spherical symmetry im- 
plies circular symmetry on the sky. The lensing angle, a(8) 
must be purely radial, dependent only on 8, the angle on the 
sky from the center of the mass distribution. Furthermore 
circular symmetry implies (e.g., NB99) that 



a(9) 



AGM{8) D LS 
c 2 8 D L D S 



(3) 



where M{8) is the mass enclosed within a circle of angular 
radius 6, and the various D's are angular diameter distances. 
The subscript L refers to the lens, i.e., the cluster, and the 
subscript S to the source, i.e, the object being lensed. For 
reference, we write down expressions for these distances as a 
function of redshift. We follow the notation of Peebles (1993) 
for cosmological quantities, except that we use the subscript 
M rather than to denote the mass density parameter. For 
Dl (Ds) we have 



* The quotation marks are meant to indicate that p c is not nec- 
essarily the density at r = 0, which may in fact be infinite, but 
rather is a characteristic density in the inner regions of the pro- 
file. Its precise meaning depends on the mathematical form of the 
profile. 




I 1 1 I 1 1 1 1 I 1 

500 1000 1500 
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2000 



Figure 1. Mass density, in units of the critical density, plot- 
ted as a function of radius for three clusters each with M200 = 
10 15 h~ 1 M(? J . One cluster has an SIS density profile, one has an 
NFW profile, and one has a KKBP profile. The cluster structural 
parameters r a and <5 C are chosen to make the three clusters as 
similar as possible. See text for further details. 



D(Z) = ±.- V{Z) 



Ho l + z 



(4) 



where z is the lens (source) redshift, and y(z) is given, in 
the case of a flat (Om + Qa ~ 1) universe, by 



y{z) 



dz' 



A) v /fi M (l + 2') 3 +°A 

The lens-source angular size distance is 



H l + z q 



(5) 



(6) 



where we have introduced the abbreviations y L = y(z L ), 

Vs = y( z s)y and Vls = Vl ~ Vs ( a g am . assuming a flat 
universe). Using the above expressions, we have 



Dls _H (1 + z l 



Dl 



DlDs 



(7) 



c VlVs 

This last form of the distance ratio will be useful in what 
follows. 

Adopting the two-parameter density profile (Eq. 1), the 
projected mass within an angle 6 on the sky is given by 



[D L 9 

M{6) = 2tt / rdr I 



where 

9 a = t s /D l 



>/fl 



^•\/r 2 + z 2 ^j dz 
ydy Tj ( K \/y 2 + x 2 ^ dx, (8) 



(9) 
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is the halo scale length in angular units. Combining Eqs. 3, 
7, and 8, and using the definition of the critical density, we 
obtain, after some algebra, the deflection angle: 



a{0) = 6 8 C 9 2 s y(z L , z s ) (3(9/9,) , 



where 



j3(x) = J ydy J r) y 1 + w 2 ^ dw 



and 



l + z r 



(10) 



(11) 



(12) 



Eqs. 10 and 11 determine the lensing angle for any adopted 
(spherically symmetric) density law. Note that the lens and 
source redshifts enter only through the multiplicative factor 
y(z L , z s ), while the shape information is contained entirely 
in the function f3(x). 



2.1 Convergence and Shear 

Our formula for a(9) enables us to calculate two other im- 
portant quantities, the convergence k and shear 7. To obtain 
the former, note that 



k — -V ■ a 



da a 

-T77 + 



(13) 



where the last equality follows from circular symmetry. Us- 
ing Eq. 10 in Eq. 13, we obtain 



K(x) = 35 c 9 3 y(z L ,z s ) [ary(z)] , 



where x = 9/9 s - 

The complex shear, 7 = 71 + yy2, is given by 

1 



71 = 2 (v>, 

and 

72 = 1p,xy 



■ i>,w) 



(14) 



(15) 



(16) 



where ip(9) is the lensing potential (NB99) defined by a = 
Vi/>, and the comma signifies partial differentiation. Circular 
symmetry implies that a = ip' , giving, after some manipu- 
lation, 

1 ..rda al 



71 = — cos 
and 



72 = - sin 2(f> 



dQ 9 

da a 
~d~9 ~ 9 



(18) 



where (f> is the angle between the radial direction and the 
x-axis. The shear amplitude is 



111 = VlT+yl = (a - a/9)/2 

= SS c 9 3 y(z L ,z s ) [x^p'ix) - 2x- 2 p{x)] , 



(19) 



where as before x = 9/9 s and Eq. 10 has once again been 
used. 

Note that Eqs. 17 and 18 imply that for (f> = (i.e., 
a point along the a-axis), 71 = j— y| and 72 = 0. However, 



because of the circular symmetry of the problem, there is 
no natural orientation for the two Cartesian axes. Thus, we 
may reinterpret the above equations as saying that at any 
point, we may adopt a local Cartesian coordinate system in 
which one axis points along the radial direction, the other 
along the tangential direction. With respect to such a sys- 
tem, the real part of 7 is always equal to |— y | as given by 
Eq. 19, while the imaginary part vanishes. As described be- 
low, we can measure complex galaxy ellipticities relative to 
such local axes. The resulting simplification in the probabil- 
ity distribution of observed ellipticities is described in § 3. 



2.2 Discussion 

The formalism described above shows that the relevant 
quantities for a lensing analysis — the bending angle a, the 
convergence k, and the shear 7 — are determined (up to cos- 
mological factors) by three quantities: the dimensionless cen- 
tral density S c , the scale radius r s , and the dimensionless 
projected mass function /3(x) given by Eq. 11. The last quan- 
tity depends only on the assumed form of the mass distribu- 
tion (SIS, NFW, KKBP, or other); the former two may then 
be viewed as parameters to be determined by the lensing 
analysis. In Appendix 1, we derive expressions for (3(x) for 
SIS, NFW, and KKBP. 

In Figure 2, we compare the resultant forms of a(9), 
k(9), 7(#), and M(9) for a fiducial cluster of virial mass 
10 1b /i _1 Mq lying at 2 = 0.5. The parameters 8 C w 27000 
and r s = 200ft -1 kpc were adopted for the NFW profile, 
while for SIS we took the core radius to be 10/i -1 kpc. The 
parameters for KKBP were determined by requiring that the 
logarithmic slope of the density profile, dlnp/dlnr, take on 
the value — 2 at r = 200/jT 1 kpc, where the NFW logarithmic 
slope has that value as well. These conditions plus fixing the 
virial mass at 10 15 'K~ 1 Mq uniquely determines the profiles 
and makes the comparison well-defined. 



3 MAXIMUM-LIKELIHOOD METHOD FOR 
DETERMINING HALO PARAMETERS 

3.1 Preliminaries 

The method is based on the (complex) ellipticities, e, of 
galaxy images, computed from the intensity moments Q a b, 

Q xx Qyy 2iQxy 



e = e x + i£ y = 



Qxx ~\~ Qy 



where 



Qab 



J I abdydx 
J Idydx 



(20) 



(21) 



and I is the measured intensity. As noted above, the formal- 
ism is considerably simplified if we transform to a coordinate 
system in which the real axis (axis 1) points along the tan- 
gential direction, and the imaginary axis (axis 2) along the 
radial direction, with respect to a radial vector joining the 
lens center and the galaxy. This transformation is accom- 
plished by a simple rotation, 
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Figure 2. Lensing angle a, convergence re, shear amplitude 7, 
and enclosed mass M plotted as functions of 8, the angle on the 
sky from the center of the mass distribution, for a cluster of virial 
mass 10 15 /i — 1 Mq at a redshift of 0.5. Each quantity is shown 
for each of the three density profiles studied in this paper: SIS, 
NFW, and KKBP. The three profiles were "matched" as described 
in the main text. Note that the parameters are the same as those 
in Figure 1. 



sin 2(f) cos 2(/> 
— cos 2cf> sin 2(f> 



(22) 



where <j) is the angle between the original z-axis and the real 
axis at the position of the galaxy. 

The ellipticity components ei and e 2 transform under 
lensing in a particularly simple way. As shown by Miralda- 
Escude (1991, hereafter ME91), the ellipticity components 
of the lensed image (subscript "L") are related to those of 
the source image (subscript "S"' ) by the transformations 

3 + esi 



and 



£L2 



1 + Sesi 



(i-O 



2\l/2 



CS2 



(23) 



(24) 



1 + Sesi 

In Eqs. 23 and 24, 5 is the distortion denned by (e.g., ME91) 



S = 1 , 2 - 

l + qi 

where 

l-(* + 7) 



9c 



1 



7) 



(25) 



(26) 



and k and 7 are the convergence and shear amplitude, re- 
spectively. 

Our method requires that we invert these transforma- 
tions to obtain the source ellipticity as a function of the 
image (lensed) ellipticity. This yields 



£S2 



-5 + em 
1 — Se L i 
{l-5 2 ) 1/2 e L2 



1 - Se L1 ' 

The Jacobian of the inverse transformation is given by 



3 



1 - 8e L1 



(27) 
(28) 

(29) 



3.2 Probability Distributions 

We seek an expression for the probability that the image 
of a background galaxy will exhibit ellipticity components 
(ei, £2) (we drop the subscript L in what follows). Following 
the work of Seitz & Schneider (1995,1996; see also Schneider 
& Seitz 1995), we write the probability density in terms of 
the image ellipticity components as 



p(ei, e 2 ; 5) = <j>[es{ei,e 2 ; S)] x J(ei; 5) . 



(30) 



where e s = y/e si t <= S2 . ■ 
the inverse transformation, Eqs. 27 and 28. The function 
4>(cs) is the distribution of intrinsic source ellipticities; it 
only depends on amplitude because of the assumed isotropy 
of source galaxy orientations, and is normalized by 



2tt 



dc ( 



(31) 



A successful application of the maximum-likelihood method 
requires an accurate model for 0(e). We discuss this further 
in § 6. 

The data set with which we maximize likelihood consists 
of TV lensed galaxy images, each with a measured complex 
ellipticity ei. We define the likelihood for the data set as 



JV 

£ 

i=l 



logp(ej; Si(S c , 9 S ), Qi) 



(32) 



where the individual probabilities are given by Eq. 30. Here 
we have explicitly indicated the dependence of the proba- 
bility on 6i, the angle between the radius vector connecting 
the cluster center to the galaxy and the a;-axis, which is 
needed for calculating ei and e 2 from the directly measured 
components e x and e y (cf. Eq. 22). The likelihood is maxi- 
mized with respect to the cluster structural parameters S c 
and 6 3 = r s /D L - 



4 SIMULATING LENSED BACKGROUND 
FIELDS 

In order to quantify the feasibility of the described method, 
we applied it to synthetic data sets, that covered a wide area 
in parameter space. These simulations are described in this 
section; the results are presented in the following section. 

The simulations involved generating 10 fields of galax- 
ies, 10' on a side with a resolution of 0.3" per pixel. Each 
field had approximately 1500 galaxies randomly distributed 
per field with an intrinsic ellipticity distribution given by 
ME91 
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0.2 0.4 0.6 0.8 1 



u< 
ft 




0.2 0.4 0.6 0.8 

q 

Figure 3. The intrinsic ellipticity distribution. The solid line 
shows p(q) (Eq. 33) , the dotted line shows the corresponding 
p(e) while the dashed line shows the observed distribution p bs( e ) 
(Eq. 35). An arbitrary normalization is chosen. The effects of pix- 
elization, increasing the number of galaxies with low ellipticities 
is also evident. 



p(q)dq oc (1 + cq) exp 



1 — exp f 








92 / 



dq. 



(33) 



where q is the ratio of minor to major axes, q = b/a and with 
c = 0.2, qi = 0.4 and q^ — 0.02. This probability distribution 
is shown in Figure 3. The choice of this functional form 
and the particular parameters for the ellipticity distribution 
is a best fit (ME91) to the ellipticity distribution of the 
deep CCD survey (Tyson 1988, Tyson & Seitzer 1988) while 
the number of galaxies was chosen such that the number of 
galaxies per unit area is characteristic of real images. Each 
of these 10 fields were assumed to be at a redshift z = 2.0 in 
a fi m = 0.3, Qa = 0.7, h — 0.7 universe. We added uniform 
Poisson photon noise to each of the images, but ignored the 
effects of a redshift distribution and seeing for simplicity. We 
realize that the latter should be included at a later stage. 
An example of a background field is shown in Figure 4. 

Each of these fields was lensed by dark matter haloes 
of virial masses 5.0, 7.5, 10.0 and 15.0 x 10 14 M Q each lo- 
cated at redshifts of z = 0.2, 0.5 and 0.8. For each halo 
of a given mass, we simulated three different density pro- 
files, SIS, NFW and KKBP, to study variations with mass, 
redshift, and profiles. The choice of the particular profiles 
enabled us to test the algorithm's discriminatory power at 
varying scales since the KKBP profile differs from NFW only 
at small radii while both differ from SIS at all length scales. 

However, since all the profiles are two parameter mod- 
els, specifying the mass does not uniquely determine the 
values of both parameters. Recognizing that the halo scale 
length (for NFW and KKBP) and the cluster core radius 



Figure 4. A sample field of background galaxies, seen as they 
would be in the complete absence of gravitational lensing. 



(SIS) are parameters that lack significant physical signifi- 
cance, we fixed those parameters in our simulations. The 
scale lengths for SIS and NFW were arbitrarily set to 10 kpc 
and 200 kpc respectively. To determine the scale length for 
KKBP, we matched the radius at which d log p/d log r = —2 
with that for NFW; this fixed the tkkbp = tnfw /\/1.8 = 
149.07 kpc. This determined the second parameter, the over- 
density 8 C or the velocity dispersion a v , for a halo of a given 
mass uniquely. 

The actual lensing used an algorithm described in Wil- 
son et al.(1996); it involved solving the lens equation 



<6) 



(34) 



where j3 is the radial position of the unlensed image, 9 is 
the radial position of the lensed image and a{6) is the lens- 
ing angle, for every point on the image frame and mapping 
it to the source frame, thereby correctly handling multiply 
imaged sources. Examples of these lensed images are shown 
in Figures 5, 6 and 7. The frames were analyzed using the 
FOCAS image software (Valdes, 1982) which extracted both 
the positions and isophotal moments required for the mass 
inversion. These were then analyzed using the method out- 
lined previously, the results of which are now discussed. 



5 RESULTS 

5.1 Implementation 

Before proceeding with the results of the simulation, we 
briefly describe the actual implementation of the maximum 
likelihood method, whose formalism was discussed in § 3. In 
what follows, we assume that the images have been suitably 
processed with extracted galaxy positions and ellipticities. 

An essential element of the maximum likelihood routine 
is the intrinsic ellipticity distribution of galaxies. Instead of 
directly using the distribution we used to generate the back- 
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Figure 5. The same field lensed with an SIS halo of mass m 
1.0 x 1O 15 M at z = 0.5. 



■ 

Figure 7. The same field lensed with an KKBP halo of mass 
m = 1.0 x 1O 15 M at z = 0.5. 



Figure 6. The same field lensed with an NFW halo of mass 
m = 1.0 x 1O 15 M at z = 0.5. 



ses) grid assuming the three (NFW, SIS, KKBP) models in 
turn. For these three types of analyses, an expectation value 
for the mass enclosed in a radius of 1.5 Mpc was computed 
by weighting the mass obtained at every point on the S c - 9 S 
grid by its corresponding probability, summing and dividing 
by a suitable normalization constant, 

{Mi50o) = — ZHM) — (36) 

where M1500 is the mass at 1.5 Mpc and p(S c ,d 3 ) is the 
probability function. We chose a fixed radius for the mass 
computation to compare the accuracy of mass estimation 
even if the incorrect model was assumed for the analysis. 
In the case where the virial mass is desired, an appropriate 
formalism is discussed in Appendix B. 

The parameter estimation and profile identification 
used more direct methods. The location and value of the 
maximum of the likelihood function determined the param- 
eters (S c ,0 s )ml- The profile type was determined by choos- 
ing the model with the greatest probability. The following 
sections summarize the results thus obtained. 



ground fields (Eq. 33), we simulated 10 background fields of 
galaxies (independent of those used in the actual simula- 
tions) and used a fit to the observed ellipticity distribution 
as the intrinsic ellipticity distribution. The best fit obtained 
was (see also Figure 3) 

P(e s ) de s oc exp(-1.75866(e s + 1.1388) 2 ) (35) 

with the constant of proportionality determined by the nor- 
malization of the probability function. Note that this is not 
directly related to Eq. 33 since pixelization and other obser- 
vational effects have been convolved into it. We discuss this 
in more detail later in this section. 

Each lensed image was analyzed by constructing the 
likelihood function (Eq. 38) on a 8 C -8 S (u v -6 3 for SIS analy- 



5.2 Profile Identification 

The profile identification results are summarized in Table 1 
and Figures 8 to 11. We found it instructive to consider the 
identification of all three profiles separately as well as the 
ability to distinguish between the SIS and the NFW/KKBP 
profiles jointly. This was motivated by the similarities be- 
tween the NFW and KKBP profiles, which differ only in the 
inner core, an area where weak lensing is not very sensitive. 

As is apparent from the data, our method appears to 
correctly distinguish SIS clusters from the NFW/KKBP 
joint profile; the correctly identified fraction being approx- 
imately 80% for clusters at z — 0.8. At low redshifts, both 
NFW and KKBP clusters are misidentified as SIS while at 
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Figure 8. The fraction of SIS clusters correctly identified as such. 
The masses are all xI0 15 M Q . 



Table I: Profile Identification Fractions 



Profile Type Redshift Mass (10 15 M o ) 

0.5 0.75 1.0 1.5 



SIS 0.2 1.0 1.0 0.8 0.9 

0.5 0.9 0.7 0.9 0.8 

0.8 0.8 0.8 0.8 0.8 



NFW 0.2 0.5 0.4 0.2 0.4 

0.5 0.6 0.3 0.3 0.5 

0.8 0.5 0.1 0.5 0.7 



KKBP 0.2 0.5 0.6 0.7 1.0 

0.5 1.0 1.0 1.0 0.7 

0.8 0.8 0.9 1.0 0.7 



NFW+ 0.2 0.5 0.5 0.55 0.85 

KKBP 0.5 0.85 0.8 0.9 0.85 

0.8 0.7 0.65 0.95 0.95 



0.5 0.75 1.0 1.5 



Mass = 0.5 



Mass = 1 .0 



HI 



Mass = 0.75 



Mass = 1 .5 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 
Redshift 

Figure 9. The fraction of NFW clusters correctly identified as 
such. The masses are all xlO 15 M0. 
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Figure 10. The fraction of KKBP clusters correctly identified as 
such. The masses are all xlO 15 Af0. 



Table 1. The fraction of dark matter profiles identified correctly, 
at the specified redshifts and masses. The NFW+KKBP category 
is the fraction of clusters correctly identified as not being SIS. 



Mass = 0.5 



Mass = 1 .0 



Mass = 0.75 



Mass = 1 . 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 
l-ie d s h i f I 

Figure 11. The fraction of clusters correctly identified as not 
being SIS. The masses are all xl0 15 M Q . 



higher redshifts, NFW clusters tend to be misidentified as 
KKBP. This implies that our method is, as expected, unable 
to distinguish between the NFW and KKBP profiles, espe- 
cially at high redshifts. This can be traced to the fact that 
the two profiles only differ in their innermost regions, evi- 
dent from Figure 1. However, weak lensing is not sensitive 
to these regions of the cluster, and so to distinguish between 
such clusters, we must fall back to strong lensing results. We 
discuss this further in § 6. For clusters at low redshifts, the 
lensing data does not extend to a large enough radius to 
be sensitive to the change in slope of the NFW and KKBP 
profiles. Accordingly, a greater number are misidentified as 
SIS. 
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Table II: Mass Estimation - I 



Profile Type 


M 1500 (1O 15 M ) 




Redshift 








0.2 


0.5 


0.8 


SIS 


0.575 


0.549 


0.550 


0.545 




0.754 


0.715 


0.727 


0.731 




0.914 


0.901 


0.875 


0.873 




1.197 


1.186 


1.191 


1.169 


NFW 


0.550 


0.515 


0.511 


0.510 




0.756 


0.710 


0.724 


0.703 




0.951 


0.906 


0.882 


0.922 




1.318 


1.298 


1.291 


1.223 


KKBP 


0.546 


0.512 


0.503 


0.502 




0.675 


0.646 


0.628 


0.625 




0.800 


0.757 


0.741 


0.742 




1.037 


0.990 


1.009 


0.957 






0.2 


0.5 


0.8 



Table 2. The mean estimated values for M1500 if we assume 
that the profile was correctly identified. 
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Figure 13. The ratio of the estimated mean mass to the true 
mass of the 10 simulated NFW clusters against the true masses. 
The errorbars are the standard deviation scaled by l/\/l0. The 
solid line is y = 1, depicting equality between the actual and 
measured masses. 
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Figure 12. The ratio of the estimated mean mass to the true 
mass of the 10 simulated SIS clusters against the true masses. 
The errorbars are the standard deviation scaled by 1/%/TO. The 
solid line is y = 1, depicting equality between the actual and 
measured masses. 




6 8 10 

M isoo (xlO 14 M s ) - Actual 

Figure 14. The ratio of the estimated mean mass to the true 
mass of the 10 simulated KKBP clusters, against the true masses. 
The errorbars are the standard deviation scaled by l/\/lO. The 
solid line is y = 1, depicting equality between the actual and 
measured masses. 



5.3 Mass Estimation - I 

The mass estimates for each of the clusters are summarized 
in Table 2 and Figures 12, 13 and 14. As mentioned earlier, 
these estimates correspond to the mass of the cluster at 1.5 
Mpc = Afisoo- Also, we assumed that the cluster profile was 
identified correctly; mass estimates for incorrectly identified 
clusters are discussed in the following section. 

Table 2 shows that the cluster masses are correct to 
within a 10% error at M1500 ~ 0.5Mq that improves to 
< 5% at M1500 « 1.5Mg. The dependence of these estimates 
on redshift appears to be weak and there is no clear trend. 
This can possibly be attributed to the fact that although 
the physical extent of the images increases with increasing 



redshift, the signal falls off with increasing radial distance 
from the center of the cluster. This loss of signal with in- 
creasing redshift is also evident from an examination of the 
diagonal elements of Table 6. 

An interesting feature of the data presented is the sys- 
tematic 5% bias in the mass estimates. Further numerical 
experiments pointed to the input ellipticity distribution as 
the cause of this bias. This dependence of the bias on the 
ellipticity distribution can be explained as follows: imagine 
we had assumed an ellipticity distribution of a 8 function at 
zero ellipticity. In such a case, all mass estimates would be 
biased above since random elliptical galaxies would be er- 
roneously interpreted as a weak lensing signal, indicating a 
mass distribution where none exists. Similarly, if we had as- 
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Table III: Mass Estimation - II - SIS clusters 



Redshift 


Mi5oo(lO 15 M ) 


Profiles 








NFW 


KKBP 


z = 0.2 


0.575 


0.300 


0.224 




0.754 


0.425 


0.311 




0.914 


0.541 


0.404 




1.197 


0.755 


0.590 


z = 0.5 


0.575 


0.382 


0.285 




0.754 


0.522 


0.405 




0.914 


0.640 


0.525 




1.197 


0.906 


0.780 


z = 0.8 


0.575 


0.378 


0.279 




0.754 


0.552 


0.410 




0.914 


0.669 


0.531 




1.197 


0.889 


0.758 






NFW 


KKBP 



Table 3. The mean estimated masses for the SIS clusters when 
analyzed assuming an incorrect profile. 

sumed an intrinsic ellipticity distribution that had a greater 
fraction of elliptical galaxies than actually existed, such a 
distribution would result in an underestimate of the lensing 
and therefore the mass distribution. 

It would appear, at first sight, that in our "experiments" 
we could avoid this by simply including our input ellipticity 
distribution, i.e. use Eq. 33 instead of Eq. 35. However, due 
to the finite size of galaxies and pixelization effects, the input 
ellipticity distribution and the observed ellipticity distribu- 
tion are not the same. Our results correspond to estimates 
made with fits to the observed distribution, a choice further 
justified by the fact that any real observations must rely on 
such observed distributions since the true intrinsic distribu- 
tion is unknown. This choice was borne out by numerical 
experiments that showed that estimates based on our Gaus- 
sian fit were more accurate than those from Eq. 33. 

5.4 Mass Estimation - II 

The mass estimations discussed in the previous section as- 
sumed that the profile of the cluster was correctly identified. 
It is equally important to quantify the bias in the mass es- 
timation if the incorrect profile is used, especially if we con- 
sider using these cluster masses for cosmology. Tables 3,4 
and 5 and Figures 15, 16, and 17 summarize the mass es- 
timation in our simulations when an incorrect profile is as- 
sumed. We also list the average bias as a function of redshift 
in Table 6. A clear trend is the reduction in the bias with 
increasing redshift; for instance the bias in KKBP clusters 
analysed as SIS clusters drops from 70% at z = 0.2 to 35% 
at z = 0.8, a reduction by a factor of two. This is easily 
explained by the fact that the degree of extrapolation de- 
creases with increasing redshift since more of the cluster is 
visible and therefore the mass is better determined. Note 
also that the bias in an NFW/KKBP family cluster anal- 
ysed with an SIS profile is of the order of 30% even at a 



Table IV:Mass Estimation - II - NFW clusters 



Redshift 


Mi5oo(lO 15 M ) 


Profiles 








SIS 


KKBP 


z = 0.2 


0.550 


0.717 


0.331 




0.756 


1.012 


0.488 




0.951 


1.253 


0.636 




1.318 


1.722 


0.991 


z = 0.5 


0.550 


0.688 


0.378 




0.756 


0.912 


0.579 




0.951 


1.120 


0.742 




1.318 


1.607 


1.110 


z = 0.8 


0.550 


0.673 


0.394 




0.756 


0.872 


0.552 




0.951 


1.111 


0.755 




1.318 


1.565 


1.096 






SIS 


KKBP 



Table 4. The mean estimated masses for the NFW clusters when 
analyzed assuming an incorrect profile. 



Table V:Mass Estimation - II - KKBP clusters 



Redshift 


M 15O o(lO 15 M ) 


Profiles 








SIS 


NFW 


z = 0.2 


0.546 


0.943 


0.734 




0.675 


1.166 


0.943 




0.800 


1.364 


1.135 




1.037 


1.738 


1.347 


z = 0.5 


0.546 


0.797 


0.611 




0.675 


0.966 


0.760 




0.800 


1.124 


0.896 




1.037 


1.433 


1.119 


z = 0.8 


0.546 


0.766 


0.582 




0.675 


0.942 


0.735 




0.800 


1.087 


0.860 




1.037 


1.331 


1.080 






SIS 


NFW 



Table 5. The mean estimated masses for the KKBP clusters 
when analyzed assuming an incorrect profile. 

redshift of 0.8, while the bias drops to 15% if the assumed 
profile was incorrect but part of the NFW/KKBP family. 
There is therefore a clear advantage to having some degree 
of profile determination before extrapolating. 

5.5 Parameter Estimation 

The parameter estimates for our simulations are summa- 
rized in Tables 7 and 8 and Figures 18 to 23. The most 
striking feature about these data is that the parameters are 
estimated to within 25% even though the masses were es- 
timated with a < 10% error. This is explained by realizing 
that there is a degeneracy in determining the parameters for 
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Table VI: Average percentage bias 



Redshift 


Actual 




Bias 








SIS 


NFW 


KKBP 


z = 0.2 


SIS 


-2.2 


-40.5 


-55.2 




NFW 


32.0 


-4.2 


-31.6 




KKBP 


70.3 


38.3 


-4.9 


z = 0.5 


SIS 


-3.1 


-28.7 


-41.6 




NFW 


19.5 


-5.2 


-20.8 




KKBP 


40.6 


11.1 


-6.1 


z = 0.8 


SIS 


-3.6 


-26.5 


-41.5 




NFW 


16.9 


-5.1 


-19.6 




KKBP 


34.9 


7.0 


-7.4 






SIS 


NFW 


KKBP 



Table 6. The observed percentage bias in M1500 estimation for 
various assumed profiles. 
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Figure 16. The estimated mean masses of the 10 simulated NFW 
clusters against the true masses, if analzyed assuming an SIS or 
KKBP profile. The errorbars are the standard deviation scaled 
by l/VTO. The solid line has slope 1 and zero intercept, depicting 
equality between the actual and measured masses. 
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Figure 15. The estimated mean masses of the 10 simulated SIS 
clusters against the actual masses, if analyzed assuming an NFW 
or KKBP profile. The errorbars are the standard deviation scaled 
by 1/vTO. The solid line has slope 1 and zero intercept, depicting 
equality between the actual and measured masses. 
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Figure 17. The estimated mean masses of the 10 simulated 
KKBP clusters against the true masses, if analzyed assuming an 
SIS or NFW profile. The errorbars are the standard deviation 
scaled by l/\/T0. The solid line has slope 1 and zero intercept, 
depicting equality between the actual and measured masses. 



a given mass, increasing the variance in the parameter es- 
timation. However, the correlated errors in the parameters 
combine to give a smaller cumulative error in the mass. 

It is also relevant to note that the errors in the esti- 
mation of the core radius for SIS clusters are significantly 
greater than the corresponding errors in the scale radius for 
NFW/KKBP clusters. This is a result of the fact that the 
core radius for the SIS clusters is ~ 10 kpc and therefore 
has little effect on the observed weak lensing. On the other 
hand, the scale radius in the NFW/KKBP clusters affects 
the mass profile of the cluster at significantly greater scales 
and is therefore probed by weak lensing. Any attempt to 
robustly estimate the core radius for SIS clusters must rely 
on strong lensing. 



6 DISCUSSION 

We have studied the reliability of mass profile determina- 
tions for clusters of galaxies using a parametric inversion 
method. We found that the mass profiles of different "fami- 
lies" of profiles could be reliably determined to within 80%. 
We appear to able to distinguish profiles motivated by N- 
body CDM simulations (in our case, NFW and KKBP) from 
the canonical r~ 2 SIS profile, making our method a promis- 
ing approach to test this prediction. However, we were un- 
able to distinguish between the various "flavours" of models 
predicted by N-body simulations since these models only dif- 
fer in their inner cores and are therefore inaccessible to weak 
lensing. Strong lensing analyses of arcs, sensitive to this in- 
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Table VII:<5 c / a v Estimation 



Profile Type 


S c (xlO 4 ) / 




Redshift 






a v (xlO 3 km/s) 


0.2 


0.5 


0.8 


SIS 


0.913 


0.886 


0.887 


0.887 




1.045 


1.016 


1.023 


1.029 




1.150 


1.140 


1.116 


1.128 




1.317 


1.310 


1.319 


1.302 


NFW 


1.57 


1.59 


1.68 


1.72 




2.15 


2.28 


1.92 


2.12 




2.70 


2.80 


3.19 


2.45 




3.75 


3.68 


3.63 


4.35 


KKBP 


2.29 


2.34 


2.41 


2.41 




2.84 


2.76 


2.93 


3.10 




3.36 


3.44 


3.50 


3.59 




4.36 


4.28 


4.03 


4.60 






0.2 


0.5 


0.8 



Table 7. The mean estimated values for S c (x 10 4 ) for NFW and 
KKBP clusters and a v (xlO 3 km/s) for SIS clusters. The profiles 
are assumed to be correctly identified. 
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Figure 18. The ratio of the estimated a v to the actual a v of 
the 10 simulated SIS clusters plotted against the true values. The 
errorbars are the standard deviation scaled by 1/%/TO. The solid 
line is y = 1 , depicting equality between actual and measured cr v . 



ner core might prove to be sufficient to constrain the inner 
parts of the profiles but this still remains an open question. 

We also examined the biases in mass estimates when 
an incorrect cluster profile (eg. SIS) is assumed and found 
that these biases can be quite significant (~ 70%), espe- 
cially at low redshifts (z ~ 0.2). At higher redshifts, the 
bias is weaker, although still substantial at ~ 50%. This 
redshift dependence was traced to the fact that at higher 
redshifts, images of a given angular extent correspond to a 
greater physical extent. Weak lensing constrains the mass 
within the observed area, but in order to estimate the virial 
mass (or Misoo), one must extrapolate to the appropriate 
radius. Using an incorrect profile for this extrapolation will 
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Figure 19. The ratio of the estimated <5 C to the actual <5 C of the 
10 simulated NFW clusters plotted against the true values. The 
errorbars are the standard deviation scaled by 1/vTO. The solid 
line is y = 1, depicting equality between actual and measured <5 C . 
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Figure 20. The ratio of the estimated 8 C to the actual <5 C of the 
10 simulated KKBP clusters plotted against the true values. The 
errorbars are the standard deviation scaled by l/\/lO. The solid 
line is y = 1, depicting equality between actual and measured <5 C . 



bias the mass estimate. However, the degree of extrapolation 
decreases with increasing redshift, reducing the bias. 

Finally, we also showed that the method could ex- 
tract the correct cluster profile parameters although only 
to within 25%. The greater errors in this case can be traced 
to a degeneracy in the parameters for a given mass; such 
a degeneracy would naturally broaden the likelihood func- 
tion and increase the errors. The case for a degeneracy is 
further strengthed by the fact that the mass estimates are 
significantly more robust, being within 10%. 

This paper represents a first attempt to demonstrate 
the feasibility of our method and estimate the reliability of 
profile determination and the bias introduced in the mass 
estimate due to an incorrect profile. For simplicity, we have 
left a number of practical aspects unaddressed which we 
now consider. The most relevant issue is the measurement 
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Table VIILS. Estimation 



Profile Type 


@s,true 




Mass 


(10 15 M Q ) 






: 


2 




f 3XCSGC.J 


U . u 


n 7^ 
u. / o 




1 F ! 

1 ,o 


= 




SIS 


4.3 


3.5 


3.6 


5.1 


5.1 


1 




1 




2.3 


2.3 


2.1 


1.5 


4.0 










1 Q 


1 Q 

1 .J 




2 4 


9 7 







8 


NFW 


86.6 


85.7 


83.5 


86.0 


97.7 


1 








46.8 


45.0 


50.2 


45.0 


54.0 


= 


: 


2 




38.0 


37.5 


38.7 


44.2 


36.2 








KKBP 


64.5 


62.7 


64.5 


62.5 


64.0 






1 




34.9 


33.5 


33.7 


33.5 


36.5 










28.3 


27.5 


26.7 


27.0 


27.0 







B 






0.5 


0.75 


1.0 


1.5 









Table 8. The mean estimated values for 8 S in arcseconds. The 
profiles are assumed to be correctly identified. 
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Figure 22. The ratio of the estimated 8 B to the actual S of the 
10 simulated NFW clusters plotted against the true values. The 
errorbars are the standard deviation scaled by 1/vTO. The solid 
line is y = 1 , depicting equality between actual and measured 8 S . 
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Figure 21. The ratio of the estimated S to the actual 8 S of 
the 10 simulated SIS clusters plotted against the true values. The 
errorbars are the standard deviation scaled by 1/vTO- The solid 
line is y = 1 , depicting equality between actual and measured 8 S . 
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Figure 23. The ratio of the estimated 8 S to the actual 8 S of the 
10 simulated KKBP clusters plotted against the true values. The 
errorbars are the standard deviation scaled by l/v / 10. The solid 
line is y = 1 , depicting equality between actual and measured 8 S . 



of the intrinsic ellipticity distribution. We computed this 
distribution self consistently in our simulation. A real world 
application must determine this background ellipticity dis- 
tribution in fields away from significant weak lensing sources. 
The Medium Deep Survey (Ratnatunga et. al 1999) is one 
such possibility. There are a number of concerns with this 
approach such as the evolution of the ellipticity distribu- 
tion with redshift and the effects of the instrument of the 
ellipticity distribution, that still remain unexplored. 

A second area that must be incorporated into the 
method is the redshift distribution of the lensed galaxies, 
since the strength of the lensing is dependent on the distance 
between the galaxy and the lens. There are two possible ap- 
proaches to solving this problem within the framework of 
our method. The first is to use photometric redshifts and 
use those as a second input parameter in Eq. 38 



JV 

C = -J2logp(e i ,z i ;S i ,e i ). (37) 

i=l 

The second is to convolve a redshift distribution directly into 
the likelihood function and measure its parameters from the 
data itself. Eq. 38 would then read 

JV 

£ = -^logp(e i ;5 i) 6l i) Zi). (38) 

!=1 

where Zi represents the parameters of the redshift distribu- 
tion. The generic effect of such an approach would be to 
broaden the likelihood function and therefore, increase the 
errors in the measurements. 

Finally, the effect of seeing on these results must be 
addressed. Seeing, as it tends to circularise images, would 
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generically reduce the lensing signal and therefore result in 
an underestimate of the mass distribution. The measured 
intrinsic ellipticity distribution would also be affected by 
seeing, and this may compensate for the effects of seeing, al- 
though a definitive answer must wait for further simulations 
that include seeing. 

In spite of as yet unexplored observational effects rel- 
evant to our technique, our formalism still calls for an ob- 
servational test of theoretical predictions of cluster profiles. 
Improvements in telescopes designs for weak lensing searches 
ought to make this a feasible approach in the future. 
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APPENDIX A: ANALYTIC FORMS OF THE 
CONVERGENCE AND SHEAR FOR SPECIFIC 
MODELS 

The convergence and shear for a circularly symmetric mass 
distribution were given by equations 14 and 19 in the main 
text . In this Appendix we derive specific expressions for the 
SIS, NFW, and KKPB cluster profiles, t 



Al SIS 

The density law for a softened isothermal sphere with a core 
radius r. is 



p(r) 



2-kGt 2 



1 + (r/r 3 y 



(Al) 



Making the connection with the formalism developed in § 2, 
r](x) = (l + x 2 )- 1 and S c = 4al/3(H r s ) 2 . This form of r](x) 
leads straightforwardly to the following expressions for (3{x) 
and its derivative: 



(3(x) 



X 2 + 1 

X 



-0 



(A2) 
(A3) 



2 VTTx 2 

and thus yields the convergence and shear for the SIS model 



k(x) 



■y(x) = 



2na 2 



-lVLS 



ys VI + x' 2 



-lVLS 



c 2 3 ys V VI +x 2 1 + VlTx 2 
where x — 9/9 3 and 9 S — t s /Dl. 



A2 NFW 

For NFW we have r](x) = x~ x (l + x)~ 2 , and thus 



(A4) 
(A5) 



P{x) 



= ydy 

Jo Jo 



dw 



\A 2 + w 2 (l + ^y 2 + ™ 2 ) 

Making the change of variable tan# = w/y, we obtain 



(A6) 



f3(x) 



ydy 
o Jo 



r/2 



cos d6 
(y + cos 9) 2 



(A7) 



We now change the order of integration, and the y-integral 
is readily done, yielding 



f 

Jo 



r/2 



cos 9 In 

f-n/2 



1 + 



cos 9 



d9 



cos 2 9 d9 
x + cos 9 



t The SIS derivation is straightforward; our derivation for NFW 
is our own, but has previously been done, albeit in slightly differ- 
ent form, by Bartclmann 1996. For KKBP, we present a numerical 
result as no analytic solution was found. 



We now write ln(l + x/ cos 9) — \n(cos9 + x) — ln(cos0), 
integrate the first term by parts, and then obtain after some 
manipulation 

/■f/2 

P{x) — Inx + 



d9 



cos 9 + x 



+ 



/•t/2 

Jo ' 



cos 9 ln(secf9) d9 - 1 



The second integral above has the value 1 — In 2 (Gradshteyn 
& Ryzhik 1994, Eq. 4.387). Thus, 

(A8) 



(3(x) = ln(x/2)+J r (x) 
where 

fTr/2 



f 

Jo 



^sTT-x' ^ 

and has the following values (Gradshteyn & Ryzhik 1994, 
Eq. 2.553): 

(i) x < 1: 



(ii) x > 1: 
1 



In 



1 - v 7 ! - ^? 



T{x) 



cos 1 (l/x) 



(A10) 

(All) 
(A12) 
(A13) 

(x + lsO) 2 ■ (AM) 

Carrying out the integral (Gradshteyn & Ryzhik 1994, Eq. 
2.553) then leads to 



v 7 ^ 2 "- 1 
(iii) x = 1: 
T{x) = 1 

Differentiating equation A8 yields 

p\x) = -+r{ X ) 



where 
T 



\x) = - f 
Jo 



r/2 



T'{X) 



Z 2 - 1 



xT(x) 

x 



(A15) 



(A16) 



and 

*-(!) = -f- 

Equations A10-A18 above enable use to calculate bending 
angle, convergence, and shear for an NFW model of specified 
8 C and 9 S , by means of equations 10, 14, and 19 of the main 
text. 



A3 KKBP 

A KKBP halo has r)(x) = x- (h2 (l + x 2 )- 1A , and thus 



(3{x) = 
where 



/ ycr 

Jo 



(y) dy , 



(A17) 
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°{V)= / (y 2 + w 2 )- (L1 [l + y 2 + w 2 ]- 1A dw. (A18) 
Jo 

Unlike the SIS and NFW halos, there is no analytical form 
for the functions /3 and a for KKBP. We have therefore 
integrated equation A18 numerically, and fit the results to 
a rational function. The following expression provides a fit 
accurate to better than 1%: 

1.250 + 0.9651331- 
a(X > ~ 1 + 1.18363a; + 1.06347a; 2 + x 3 ' { ' 

Substitution of this expression into equation A17 allows one 
to obtain an analytic expression for (3(x) via integration: 

0(x) = -0.0558209 + /3i (x) + f3 2 (x) + /3 3 (x) , ( A20) 

where: 

Pi(x) = 0.126086 tan -1 (0.0605662 + 0.970072a;) (A21) 

fo{x) = -0.176946 ln(0.937601 + x) (A22) 

(h{x) = 0.571039 ln(1.06655 + 0.124869a; + x 2 ) . (A23) 



msis(y) = V - tan 1 y (B5) 
and 

m NFW (y) =ln(l + y)- . (B6) 

It is useful to express the virial mass in terms of typical 
halo parameters. From equation Bl one finds: 

My = 3.49 x 10" 3 (A) m(yv) h -i M0 . (B7) 

For example, a KKBP halo with for r 3 = 100 kpc and S c = 
10 5 has M v = 8.7 x 10 14 ft _1 M s , and r v = 1.55 Mpc, for 
S v = 200. 



APPENDIX B: DETERMINATION OF 
CLUSTER VIRIAL MASSES 

The virial mass of a cluster, Mv , is defined to be the mass 
within a sphere whose mean interior density is larger than 
the background density by a particular factor. The radius of 
the sphere is called the virial radius, Ry- In an Einstein-de 
Sitter universe the factor is 178; in open and flat low-density 
universes the factor is larger. For the purposes of this paper, 
we simply assume the factor is known, and denote it 5v- (In 
the main body of the paper we take Sv = 200.) In this 
Appendix we show how to calculate My and Rv for a given 
density law. 

As usual, we write the halo density profile as p(r) = 
Pcrit5 c ri(x), where x = r/r 3 and r s is the characteristic ra- 
dius. The mass interior to a radius R is thus 



M(R) = 4nrsp C rit5 c m(y) , 



(Bl) 



where y — R/r s and m(y) is the dimensionless mass profile 
given by 



m(y) 



x r)(x)dx ■ 



(B2) 



The condition that the mean density within Rv is Sv 
times the background density is 

M(r v ) 



4nr v /3 



Sv^M per it • 



(B3) 



(B4) 



Combining this with equation Bl yields 
m(y v ) = ^-n M y v , 

where yv = Rv/r 3 . For any given density law p(r) and 
value of Q.m, equation B4 can be solved for yv, which in 
turn yields the virial radius Rv = yvr a , and the virial mass 
through equation Bl. 

The integral in equation B2 is readily carried out for 
the SIS and NFW density profiles, yielding 
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